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The ability for most hadrons to decay via strong interactions prevents the direct measurement of 
their electromagnetic properties. However, a detailed understanding of how these resonant states 
feature in scattering processes can allow one to disentangle such information from photo production 
processes. In particular, there has been increasing interest in the determination of magnetic dipole 
moments using such methods. In a recent study [J], Gudino et al. provide the first experimental 
determination of the magnetic dipole moment of the rho meson. To facilitate a comparison with this 
experimental determination, we present a calculation of the rho meson and pion electromagnetic 
form factors calculated in the framework of Lattice QCD. Using the PACS-CS 2+1 flavour full 
QCD gauge field configurations, we are able to access low values at near-physical quark masses. 

Through the use of variational techniques, we control excited state systematics in the matrix elements 
of the lowest-lying states and gain access to the matrix elements of the first excited state. Our 
determination of the rho meson g-factor gp = 2.21(8) is in excellent agreement with this experimental 
determination, but with a significantly smaller uncertainty. 


I. INTRODUCTION 

Understanding how QCD gives rise to the rich and di¬ 
verse structure of the hadrons remains an ongoing effort 
both theoretically and experimentally. Probing the elec¬ 
tromagnetic structure remains one of the most effective 
methods. By measuring the electromagnetic form factors 
of these states, we map out the distribution of charge and 
magnetism within and gain valuable insight into the un¬ 
derlying dynamics. 

Unlike the nucleon and pion, the vast majority of QCD 
eigenstates are unstable with respect to the strong in¬ 
teraction. This makes experimental determinations of 
their properties via conventional means impractical due 
to their short lifetime. Instead one is forced to disen¬ 
tangle resonant state properties from radiative processes 
in scattering experiments. Recently, there has been in¬ 
terest in extracting magnetic dipole moments [IHi. In 
particular, Gudino et al. [l| have been able to provide 
a determination of the magnetic dipole moment of the 
rho meson, the first such experimental result for a vector 
meson. 

It is thus timely to consider what results lattice QCD 
can provide as a comparison. The overwhelming focus 
within the lattice community has understandably been 
on the nucleon, and to a lesser extent the pion, where 
there exists a large body of experimental data [^. In 
comparison, there have been oi^ a handful of lattice 
studies of rho meson structure [5l-l7| with the majority 
of these being performed in the quenched approximation 
to QCD. The only existing calculation with dynamical 
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quarks was a preliminary calculation Q by the QCD-SF 
collaboration using moderate to heavy pion masses and 
thus insensitive to the chiral dynamics of full QCD. It is 
also limited to a large value of nontrivial 4-momentum 
transfer with = 0.44 GeV^. Thus, there is an ur¬ 

gent need for a determination of the rho meson form fac¬ 
tors at low and near-physical quark masses in full 
QCD to allow for a direct comparison with experiment. 

The progress of the past few years had lead to lattice 
studies near or at physical masses and on lattice volumes 
significantly larger than the states that occupy them. 
Working in this new regime has reduced the need for 
chiral extrapolation, but there are new systematics that 
must be considered to ensure accurate results. As quark 
masses decrease, signal weakens and so one is forced to 
work in early Euclidean-time regions where excited state 
contamination can be significant. Several studies have 
explored new approaches to properly handle excited state 
contaminations with the most popular of these being the 
summation method ii and the variational approach 
[IMll- Here we make use of the variational approach. 

In our recent work [T^, we demonstrated how this 
method provides improved plateau quality and early on¬ 
set of ground state dominance without the need for the 
fine-tuning of source and sink operators. Furthermore, 
this improvement allows us to probe the system at ear¬ 
lier Euclidean times giving rise to significant improve¬ 
ment in the statistical uncertainties. Preliminary studies 
of electromagnetic properties provide similar con¬ 

clusions. 

In this work we present a calculation of the rho-meson 
and pion form factors utilising variational methods al¬ 
lowing for an accurate determination of these quantities 
at near-physical masses. Furthermore the variational 
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approach gives us access to the form factors of excited 
states. The use of moderately large lattice volumes al¬ 
lows us to extract our results at low 

For the extraction of the form factors we follow the 
electromagnetic form factor formalism of Hedditch et al. 
0 with the necessary changes required for use within the 
variational approach. 

The remainder of this paper is organized as follows. In 
Section 2 we discuss the use of correlation matrix tech¬ 
niques in the evaluation of hadron matrix elements. In 
Section 3 we outline the framework used to extract the 
pion (pseudoscalar) and rho-meson (vector) form factors 
on the lattice, while Section 4 summarises the details of 
our lattice simulation. In Section 5 we present our results 
and finally in Section 6 we provide concluding remarks. 

II. THEORETICAL EORMALISM 

On the lattice, the determination of hadronic observ¬ 
ables begins with the evaluation of the Euclidean-time 
two and three-point correlation functions. Given an in¬ 
terpolator x(a;) for the state in question, the two-point 
function is defined as 

G{p,t) = ^ (fl| x(x) xt(o) |0), (1) 

X 

and the three-point function 

Go{p',pMM) = ^ 

X2,Xi 

X (O|x(a;2)O(a;i)x'''(0) \il ), 

( 2 ) 

where O is the current used to probe the state. Here we 
are interested in the electromagnetic current O = J^. In 
general, this interpolator xi^) have overlap with all 
eigenstates consistent with the given quantum numbers. 
Inserting a complete set of states between our interpola¬ 
tors in Eq. © 

G{p,t) = ^e"-®^“®*(H|x(0) |a,p)(a,p|x'''(0) |H), 

a 

we can see that the contribution to the correlation func¬ 
tion from each state, a, is exponentially suppressed by 
its energy. The standard approach for examining the 
ground state is to work at large Euclidean time where ex¬ 
cited states are suppressed. In the past this approach has 
been sufficient for most quantities considered, but with 
ensembles at near physical masses and studies seeking 
precision determinations of hadronic observables, there 
has been increased concern as to whether Euclidean time 
evolution is a sufficient measure for eliminating excited 
state effects [iMl- In this study we select the varia¬ 
tional approach to isolate the individual contribution of 
a particular eigenstate to the two and three-point corre¬ 
lators. 


The variational method [2^ [2l| has become a standard 
approach for the study of the excited state spectrum of 
QCD. The goal of this approach is to produce a set of op¬ 
erators (/)“ that couple directly to individual QCD energy 
eigenstates 

(D|<^“|Ap)cxd“^ (3) 

The way in which this is achieved is to take an existing 
basis of operators {xi{x)\i = 1, ...,n} and construct 
the optimised operators (j)°‘ as linear combinations 

=^vfxi{x), (/)“^'(x) = . (4) 

i 3 

Beginning with the matrix of two-point correlation func¬ 
tions 

Gy (P,0 = E (^1 xl(0) 1^2), 

X 

and making use of Eqs. m and 0, both Gij {p, t) and 
Gij (j), t) uj have a recurrence relation governed by the 
energy of state a. It follows that the necessary vectors 
wf and are the solutions of the generalised eigenvalue 
equations 

vT Gij{p,tQ + 5t) , (5a) 

Gij{p,to + 6t)uJ Gij{p,tQ)uj. (5b) 

It is important to note that Eqs. (153 and (I5bll are evalu¬ 
ated for a given 3-momentum p and so the corresponding 
operators satisfy Eq. ([3]) for this momentum only. The 
correlator for the state \a,p) can be obtained from the 
correlation matrix by projecting with the relevant eigen¬ 
vectors 

G(p, t; a) = Vi {p) Gy (p, t) {p). 

The extension to three-point correlators follows sim¬ 
ply. Beginning with the matrix of three-point correlation 
functions 

iGohip',P,t2,h) = 

X2-,X\ 

X {^\Xii.X2)0{xi)x]{^) 1^2) 

and having evaluated the eigenvectors with the two-point 
correlation matrix, it is a simple matter of projecting 
the necessary eigenvectors onto the matrix of three-point 
functions, where care is taken to to ensure that the pro¬ 
jection is done with the correct momenta for source and 
sink, 

Go{p',P,t 2 ,ti;a) = vf{p') (Go)y (p',p, ^ 2 , G) m“(^ ■ 

o 

This applies equally well to hadron transitions (L —>■ R) 
where one simply projects with relevant eigenvectors for 
the differing left (L) and right (i?) eigenstates. 
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The application of this approach to states with explicit 
spin-degrees-of-freedom is straight forward. For a state 
with spin-1, relevant to this work, the matrix of two-point 
correlation functions is 

t) = ^ x^A^) xlAo) \n) , 

X 

where the Roman indices denote the operator in the vari¬ 
ational basis and the Greek indices denote the Lorentz 
indices. Following the arguments set out above, we can 


construct eigenvalue equations for each Lorentz compo¬ 
nent and define the corresponding optimised operators 
as 

C(a;) = XiA^), At > (6) 

* j 

i.e. we determine the eigenvectors for each Lorentz com¬ 
ponent separately such that the operator maximally iso¬ 
lates the spectrum observed in each Lorentz component 
[2^. Thus our projected two and three-point correlation 
functions are 


Gar{p,t; a) = vAip) G^j^ar{p,t) U^ip) 

{Go)aT{p\p,t2,ti;a) = vAiA) iGo)ij,ar{p',At2,ti) U^A) ■ 

Having obtained the projected two and three-point correlation functions, determination of matrix elements follows 
in the standard way with the construction of a suitable ratio to isolate the desired terms. We use the ratio defined 
by Hedditch, et al. 


RArA'A'^a) = 


' {GAA'A,t2,ti;a)) {GAAA',t2,ti;a)) 
{G^^A'(GrrA^Aa)) 


( 7 ) 


where repeated indices are not summed over. We note 
that the use of the square root requires the sign of the 
result to be recovered from the individual three-point 
functions. Moreover the covariant/contravariant place¬ 
ment of indices is for clarity only. The advantage of this 
construction is that we have exact cancellation of the 
momentum-dependent Z°‘{p} factors. In the case of the 
pseudoscalar pion interpolating field, the err indices are 
not required and can be ignored. 

As in Ref. d, [H, HI] we consider both {U} and {U*} 
configurations in creating an improved unbiased estima¬ 
tor |25| . This requires the evaluation of both +q and 
—q SST-propagators, thus incorporating parity symme¬ 
try. A consequence of this is that our correlators are 
perfectly real or imaginary, depending on the matrix el¬ 
ement under consideration. As established in Ref. [2^, 
this approach provides improved plateaus and a reduc¬ 
tion of statistical uncertainties beyond a naive doubling 
of the number of configurations considered. 


III. MESON FORM FACTORS 
A. Pseudoscalar Mesons 

Having defined suitable operators for creating and iso¬ 
lating a particular state a, it now stands to understand 
how the form factors are extracted from the ratio defined 
in Eq. ©• We begin this discussion with tt meson. 

As pseudoscalar mesons are spinless particles, the op¬ 


erator overlap can be decomposed simply 

(H| A’^{o) \TTf,{Ai)= A , 

A2EaA 

where 2“^ A) enumerates the coupling strength of the op¬ 
erator A to the state | (p)). It follows that the two- 

point function can be expressed as 

p-Ea{p)t2 

GA, t 2 ; a) = ZAp') Z‘^Hp) ( 8 ) 

and the three-point correlation function 


G'"A'A>h,ti;a) = 


-Ec<,p') -Eo,{p)ti 


2AeaA)eaA 

xzAp')z‘^\A ( ^Ap) I I MA ), 

( 9 ) 


with the matrix element {'KaA')\J^A)\'EaA)) encoding 
the interaction vertex of the pion with the electromag¬ 
netic current. For a pseudoscalar meson, this vertex is 
parametrized by a single form factor Gc{Q^) 


AAp') I J^A)\'EaA) = 

1 


2AEAp')EAA 


[A + P^]GAQ^), (10) 


where p and p' are the incoming and outgoing 4-momenta 
and = —q^ is the space-like momentum transfer. Sub- 
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stituting Eqs. ([8|), (H]) and (fTOl) into Eq.©, we obtain 

1 


which satisfies the spin sum relation 


R^[p',P]ot) = 


2^E^{p')E^{p) 


[p^+p'^]GUQ^). 


Together with our choice of kinematics, p' = {Ea,Px, 0,0) 
and p = (toq,, 0,0,0) with E^ = \/+ p ^, we access 
Gq{Q‘^) through 


TO _ 2-\/?71 q,£/q, g 


GUQ') = 


Ea + ma 


R^iPx,0;a). (11) 


J2e:{p,s)erip,s) 


= -\ 9. 


PaPr 


( 12 ) 


B. Vector Mesons 

For a vector meson, the overlap term now includes the 
spin polarization vector e(p, s) 

cafi Using these expressions, the two-point function takes the 

(111 C’^(O) I PpiP, s )) = . 2:“(^ e“(p, s), form 

,y2Ea{p) 




(13) 


while the three-point correlation function is 


G>^^{p',p,t2,ti;a) = 




Z“(p)Z“t(p)e;“(p',s') (p„(p',s')| J'^(0)|p„(p,s)) er(PiS), (14) 


where again the matrix element (pq,(p s')| J'^(0 )|pq,(p, s)) contains all of the information governing the interaction 
with the electromagnetic current. For a spin-1 system, the interaction vertex can be described by three independent 
vertex functions l2bl-l^ 


{p^{p',s')\J^{Q)\pMs)) = 


1 


2,/KWKW) 


where 


er{p',s') T^^^p\p) e“(p,s) 


r^^^{p',p) = - I 5'"^ K +p'n Gi(Q2) + G2{Q^) - G^iQ^) 

Taking the appropriate linear combination of these vertex functions, we arrive at the Sachs decomposition 

Gq(Q2) = Gi(Q2) _ G2(Q2) + (1 + 7?) GaiQ ^), 

Gm(Q') = G2(Q2), 

Gc(Q') = Gi(Q2) + ^^Gq(Q2), 

q 2 

where p = 

Substituting the vertex into Eq. (fT4l) and applying Eq. m, the three-point function becomes 

p-Ea(p')(t2-H) „-Ec,(p)tp 

G'^r(.P'^Et2,ti;a) = - —===== - Z°‘{p) Z°‘\p) A„^r{p' :P) , 

^V^AP )Ra{p} 

where we have grouped all covariant indices into the term 


\ / \ / 


(15) 
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Using the symmetry property Aa^ t{p'^ p) = AAa{p,p') and our choice of kinematics, the ratio, defined by Eq. (O, 
becomes 


Ra^T{Px,0) = 


P'.Po 

2 y/maEa V 


-1/2 


Aa'"r{Px,0)- 


In particular we consider 


Ai\ip,,0) = ^(iE^+m^)GciQ^) + l^GQiQ^)) , 
rrio, \ o rua ) 

^ 2 ° 2 (p., 0) = ^3°3 (p., 0) = f (U„ + m„) Gc(Q') - Gq(Q2)') , 

V 3 ma J 

A'3(p.,0) = -—p.Gm(Q'), 

nia 

A3^l{Px^ 0 ) = -\-px Gm{Q^) ■ 

As such, the form factors for the rho meson are isolated from the following combination of ratio terms 

(16) 

(17) 

(18) 


Px 

, = ^2R,\ - - i?3°3) . 

pI 


C. Extracting Static Quantities 

The Sachs form factors Gc, Gm and Gq describe the 
distribution of charge, magnetism and charge asymme¬ 
try within the hadron. In particular, the value of these 
functions at = 0 define the hadron’s total charge, q, 
magnetic moment, p, and quadrupole moment, Q, 

q= e Gc(0 ), (19a) 

P=:^Gm(0), (19b) 

2m 

Q = ^Gq(0), (19c) 

where m is the mass of the hadron. In this calculation, 
we work with fixed-current SST-propagators. This allows 
us to explore the form factors of many different hadrons 
without the need for additional inversions, but limits us 
to a single 3-momentum transfer and thus a single for 
each quark mass. 

We make use of a monopole ansatz for the depen¬ 

dence 

= (20) 

as suggested by a Vector Dominance Model hypothesis. 
Here we use a conserved current, Gc(Q^ = 0) = 1, and 
so we extract the monopole squared-mass via 


Motivated by the observed scaling behaviour for Ge and 
Gm of the nucleon at low , we shall assume the meson 
sector displays similar scaling for each quark sector and 
use this to extract a value for Gm{Q^ = 0) 

= = 

to facilitate a comparison with the experimental pre¬ 
diction of [l| and model expectations. Drawing on the 
monopole ansatz, we shall also use this for Gq. 

For the mean squared charge radius, we use the stan¬ 
dard definition from the small expansion of the 
Fourier transform of the charge distribution which gives 

(’■“) = . (23) 

Using our monopole form we have 

IV. SIMULATION DETAILS 

Two and three-point correlation functions are evalu¬ 
ated following the prescription outlined in [^.[^. [^1^ . 
For this calculation we make use of the PACS-CS (2-f-l)- 
flavour dynamical-QCD gauge field configurations 

















6 


TABLE I. Ensemble parameters used in this calculation. 


^ud 

rrir (MeV) 

F cfgs 

F grcs 

0.13700 

702 

350 

350 

0.13727 

570 

350 

350 

0.13754 

411 

350 

700 

0.13700 

296 

350 

700 

0.13781 

156 

197 

~1600 


made available through the ILDG [s^. These configu¬ 
rations are generated using a non-perturbatively 0(a)- 
iinproved Wilson fermion action and Iwasaki gauge ac¬ 
tion on a 32^ x 64 lattice with periodic spatial boundary 
conditions. The scale is determined via the static quark 
potential. The value (3 = 1.90 gives a lattice spacing 
a = 0.0907(13) fm resulting in a physical volume of spa- 
cial extent i = 2.9 fm. We have access to 5 light-quark 
masses, with the strange quark mass held fixed. The 
resulting pion masses range from 702 MeV to 156 MeV. 

A fixed boundary condition is applied to the temporal 
direction and our quark sources are inserted at tsrc = 16. 
We have verified that reflections associated with the fixed 
boundary are negligible at Euclidean times greater than 
16 time slices away from the boundary. Due to the lim¬ 
ited number of configurations at the lightest mass we 
make use of multiple quark sources on each configura¬ 
tion. This is done by using two maximally separated 
spacial sources with a relative temporal boundary shift 
of 8 time slices. The temporal boundary is then shifted 
by multiples of 16 time slices for each of these spacial 
sources. Multiple sources are also used for the next two 
lightest quark masses with a single spacial position, sep¬ 
arated temporally by 32 time slices. Table |T] provides a 
summary of the simulation details for each K-value and 
the corresponding value for m^. 

The three-point correlation functions are evaluated us¬ 
ing a sequential source technique with the current held 
fixed as outlined in [H, [23, [13 ■ For the current we 

use an (!l(a)-improved conserved vector current obtained 
using the Noether procedure with the improvement term 
constructed in the form of a total four-divergence [13, H^. 
The current is inserted at tj, = 21 relative to the quark 
insertion time at tsrc = 16. The resulting source-current 
separation would in general be too small and so give rise 
to excited state contamination, but as was highlighted 
in [H our use of the variational approach gives rise to 
rapid ground state dominance after the source allow¬ 
ing for considerably early current insertion and smaller 
source-current-sink separations. 

Given our choice of correlation-function ratio, we re¬ 
quire SST-propagators for both +q and —q where we 
choose q = ^x. The current is polarized with /r = 3,4 
as required by our evaluations of the form factors. As 
outlined in Ref. [1^, the invariance of the lattice action 
under U ^ U* implies the link variables {t/} and {U*} 
are of equal weight and so we choose to account for both 


sets of configurations in the evaluation of our correlation 
functions. Rather than performing further matrix inver¬ 
sions on the {U*} configurations, we make use of the 
fermion matrix property 

M{{U*}) = (C M{{U})C-^y 

with C = to give us the propagators for the {U*} 
from the existing {U} propagators 

5(a;,0;C/*) = {C S{x,0;U) . 

For the fixed-current SST-propagators this identity takes 
the form 

S{x,0;t,q,^i;U*) = {C S{x,0-t,-q, ^i;U) C~^)* . 

Thus, both q and —q momentum insertions are required 
as discussed at the end of Section 2. Our error analysis 
is performed using a second-order jackknife where the 
X^/dof for our fits is obtained through the covariance 
matrix. 

As the primary goal of this work is to cleanly extract 
matrix elements of the ground state, we choose to work 
with a small variational basis. Our operators are local 
meson operators of varying widths. This is achieved by 
applying increasing levels of gauge invariant Gaussian 
smearing to our quark sources and sinks. The smear¬ 
ing procedure is: 

tlji{x,t) = '^ F{x,x')ilJi-i{x',t) , (25) 

x' 

where 

F{x, x') = {1- a) (26) 

3 

"g ^ ^ uj^{x — jl) 5x',x—fi\ 1 

where the parameter a = 0.7 is used in our calculation. 
We use the four different levels of smearing examined in 
Ref. im to provide an optimal basis for these ensem¬ 
bles. Table |ll] lists the corresponding RMS radii for the 
quark sources. We focus on a single spin-flavour con¬ 
struction for the meson interpolators and draw on the 
various source/sink smearing widths to enable the effi¬ 
cient and accurate isolation of states in the variational 
approach. The use of a variety of widths enables the re¬ 
sulting operators to form nodal structures in the radial 
wave functions of the excited states [s^l and is central 
to our ability to rapidly isolate states in our two- and 
three-point correlation functions. 

For the spin-flavour form of our local interpolator we 
choose to use, 

X7r{x) = d{x)x5u{x) 

Xp{x) =d{x)xiu{x) 
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for the TT and p meson respectively. When coupled with 
our four smearing widths this allows for the construc¬ 
tion of up to a 4 X 4 correlation matrix. The use of the 
alternate bi-linear forms 7075 and 707 ^ was considered, 
but these were found to not provide any additional basis 
span when used with more than two smearing levels. We 
considered all combinations of variational parameters to 
and St in range 17-20 and 1-4 respectively where a super¬ 
position of states can be used to constrain the analysis. 
With regard to state isolation and the stability of the 
analysis, the optimal choice is to = 17 and St = 3 for 
the three heavier masses and to = 17 and St = 2 for the 
two remaining lighter masses. The use of an earlier to 
value relative to baryon studies [H, [13, on the same 
ensembles is unsurprising given the larger energy gaps 
displayed between the ground state and first excitation 
in the meson sector. 

In performing our variational analysis we choose to 
use the symmetrisation procedure outlined in Ref. [39l| . 
Namely we exploit the ensemble average symmetry 

and consider the improved unbiased estimator 
2 0 “f Gji{p^ t)] 

Enforcing this symmetry allows us to re-express the 
eigenvalue equations (I5al) and HEI in terms of a real sym¬ 
metric matrix 


[G-^/\pM)G{,pM+5t)G-^'\pMh < = . 

The resulting eigenvectors from this formulation form an 
orthogonal basis. These orthogonal eigenvectors, wf{p), 
can also be obtained from the eigenvectors ufip) using 

wf (p, to) (jP). 

Upon normalising this basis, we are able to sort the eigen¬ 
vectors consistently across jackknife complement sets, 
and track eigenstates across momenta using the approx¬ 
imate orthogonality condition 


Correlators are normalized to be 

Gij {p-jt) 


0 ( 1 ) by considering 
1 


Vg ii (P, tsrc) 


'S/Gjj (p, tsrc) 


TABLE II. The RMS radii for the various levels of smearing 
considered in this work. 


Sweeps of smearing 

RMS radius (fm) 

16 

0.216 

35 

0.319 

100 

0.539 

200 

0.778 



FIG. 1. Low-lying eigenstates in the pion (0'^ ) channel from 

our 4x4 correlation matrix of smeared sources and sinks. The 
left-most points at the physical pion mass are experimental 
measurements of the spectrum. 


The resulting eigenvectors then give us a direct measure 
of the relative contribution from each interpolator. Ac¬ 
cordingly, we also normalise the three-point correlators 
using the relevant two-point functions 


"s/Gii{j)\ tgj-c) 


Gf,(/,p,t 2 ,U) 


yGjjipptg^^ 


V. RESULTS 

A. Low lying meson spectrum 

In Figs. 1 and 2 we display the resulting spectrum be¬ 
low 3 GeV in the 0 ^ and 1 channels respectively 
obtained from the 4x4 correlation matrix of smeared 
sources and sinks. The masses for the states considered 
in the form factors analysis are summarised in Table IIIII 

Beginning with the pion channel we find three well 
separated eigenstates, consistent with the spectrum ob¬ 
served in previous works examining the entire isovector 
meson sector [40l - l43| . At the lightest mass we find that 
our first excited state is consistent with the 7r(1300). A 
notable feature in our spectrum is a significant shift in 
the mass for the first excited state at the second light¬ 
est mass and a similar jump for the second excited state 
at the middle mass. This feature is observed across the 
range of variational parameters considered and on cor¬ 
responding 3x3 analyses formed from subsets of the 
variational basis. A similar feature is observed in posi¬ 
tive parity spectrum of the nucleon (^ . . Examina¬ 

tion of the wave functions for these nucleon excitations 
[s^ show significant finite volume effects for the light¬ 
est two masses which may give rise to an increase in the 
eigenstate energies. It is possible that we are observing 
a similar effect here. 

In the rho meson channel, we again observe three well 
separated eigenstates. However, in this channel we ex- 

















TABLE III. Masses for the first two tt and p meson eigenstates (ground state and first excitation) for various values of the 
hopping parameter, k, extracted from our 4x4 correlation matrix of smeared sources and sinks. 


K 

ml (GeVQ 

m^ (GeV) 

mp (GeV) 

(GeV) 

mp* (GeV) 

0.13700 

0.3876(11) 

0.6226(9) 

0.981(5) 

1.66(7) 

1.68(7) 

0.13727 

0.2647(9) 

0.5145(9) 

0.917(6) 

1.59(4) 

1.72(5) 

0.13754 

0.1509(7) 

0.3884(9) 

0.867(6) 

1.55(8) 

1.77(5) 

0.13770 

0.0811(6) 

0.2848(11) 

0.832(10) 

1.77(4) 

1.81(7) 

0.13781 

0.0260(10) 

0.1613(31) 

0.793(14) 

1.50(23) 

1.77(11) 



FIG. 2. Low-lying eigenstates in the rho meson (1 ) channel 

from our 4x4 correlation matrix of smeared sources and 
sinks. The left-most points at the physical pion mass are 
experimental measurements of the p-meson spectrum. 


pect to see two eigenstates near 1600 MeV separated by 
about 250 MeV. This suggests that a basis of local op¬ 
erators is not able to isolate both of these eigenstates. 
A similar conclusion was found in Ref. [4l|. In Ref. 
they found that the p(1450) required a basis that con¬ 
tained displaced, derivative operators. Quark model ex¬ 
pectations predict that for these two states, one is S'-wave 
dominant while the other D-wa.ve dominant Q. Given 
the radial symmetry of our operators, it is not possible 
to form H-wave states and so it is not surprising that we 
are unable to isolate both states. 

The close agreement of our first excited state with the 
p(I700) and the absence of the p(1450) from our spec¬ 
trum suggests that the p(1700) is the S'-wave state while 
the p(1450) is the D-wave state. 

Dudek et al. were able to successfully isolate both 
of these states at larger quark masses and found a trend 
consistent with our identification of states. They found 
the S-wave dominated state to be the lighter state at 
their heavier quark masses, and observed that with de¬ 
creasing quark mass the mass splitting between these two 
states became smaller, with the states nearly degenerate 
at their lightest mass of m,r — 400 MeV. A continua¬ 
tion of this trend to lighter quark masses places H-wave 
dominated state lower in mass. 


An important feature in the QCD spectrum is the 
possibility of multi-particle intermediate states. In the 
infinite-volume limit this renders the majority of hadrons 
unstable under the strong interaction. However on the 
finite-volume lattice, the QCD eigenstates are stable and 
are composed of admixtures of both single-particle and 
multi-particle states. 

Some insight into the composition of states can be 
taken from the physical spectrum and scattering thresh¬ 
olds. However, the position of these thresholds change 
on the finite volume. Multi-particle states are forced to 
overlap in the finite volume, giving rise to a volume- 
dependent interaction energy. Mixing with single¬ 
particle dominated states further distorts the spectrum to 
the point where intuition from infinite-volume scattering 
thresholds and the physical spectrum becomes irrelevant, 
particularly in volumes with lattice length L ~ 3 fm. 

Below the finite-volume modified two-particle scatter¬ 
ing threshold, states are generally single-particle dom¬ 
inated but still contain important contributions from 
nearby scattering channels. The position of states in the 
spectrum can be changed by varying the quark mass or 
the volume of the lattice and the eigenstates can become 
maximally mixed making their traditional identification 
as scattering states or resonant states impossible. In the 
case where a low-lying finite-volume scattering threshold 
sits well below the resonant state, then the lowest-lying 
state may be regarded as a two- or multi-particle scatter¬ 
ing state and the single-particle dominated state is now 
the higher eigenstate. 

In the case of the vector meson under study here, we 
must be careful to ensure that the state we are excit¬ 
ing on the lattice is in fact the single-particle dominated 
resonant state. Though there is strong evidence to sug¬ 
gest that local meson operators couple poorly to scatter¬ 
ing states, especially on larger volumes such as that un¬ 
der investigation here, we perform a check to determine 
whether the eigenstate isolated in our correlation-matrix 
analysis is single particle in nature. 

For all our ensembles, the ground state p meson at rest 
is well below the tttt threshold, and will be single-particle 
dominated. However, upon applying the boost to mo¬ 
mentum g, the extracted energy eigenstate now sits above 
the lowest-lying bare tttt energy allowed by momentum 
conservation. In order to determine whether the state we 
have isolated in the boosted case is the finite-volume p- 
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FIG. 3. An example of our comparison between the disper¬ 
sion relation E = \/m? -|- np\pmin\'^ (blue line), and the en¬ 
ergies extracted from the Hnite-momentum correlators (blue 
circles) for the rho meson over a range of momenta. Here 
m is taken from the zero-momentum correlator, \'Pmin \ is the 
magnitude of the lowest non-trivial momentum on the lattice. 
The red diamonds provide the corresponding non-interacting 
TTTT-energies allowed by momentum conservation. Results at 
the other values of m-n are similar. 


meson or the lower-lying tttt scattering state, we compare 
the extracted eigenstate energy against the expected en¬ 
ergy given by the dispersion relation for a single particle. 

In Fig. [3 corresponding to the second-lightest quark 
mass considered, we overlay the dispersion relation ex¬ 
pectation (blue band) with the energies extracted from 
the finite-momentum correlators (blue data points) for a 
range of momenta. Here we can see that the energies 
extracted from the boosted correlator are in excellent 
agreement with dispersion. This is observed across all 
the masses considered. 

As a further check, we compare the non-interacting 
TTTT energies (red data points) allowed by momentum and 
parity conservation. For the single unit of lattice momen¬ 
tum relevant to our form factor analysis, we find that 
the mass separation between the dispersion result and 
the non-interacting tttt energy is significant. Moreover, 
the attractive hnite-volume interaction in the tttt system 
would act to further increase the separation between the 
single and multi-particle dominated states. As a matter 
of principle, we do expect a tttt scattering state to reveal 
itself in the long Euclidean-time tail of our correlation 
function. However, our interpolating fields have rendered 
this contribution to be negligible at the finite-Euclidean 
times considered. 

This indicates that our correlation functions are indeed 
dominated by the resonant-like state of interest and not 
the lower-lying finite volume scattering state. Similar 
separations are found for all other masses with the ex¬ 
ception of the middle mass. Through this process we 
have determined that the state we have isolated in the 
boosted system is in fact the state most closely related 


to the resonant p meson. 


B. Form Factor Determination 

A key finding from our calculation of gA using the vari¬ 
ational approach [T^, was that the optimised operators 
obtained from the correlation matrix analysis composed 
of various smeared source/sink operators gave rise to 
significant improvement in the quality and duration of 
plateaus from which the matrix elements were extracted. 
Comparison with a modest historically-typical choice of 
smearing highlighted that excited state effects are signif¬ 
icant and suppress the value of gA- 

We are able to draw similar conclusions here. In 
Fig. 0] we present a comparison of the correlation func¬ 
tion ratios providing the rho-meson charge, magnetic and 
quadrupole form factors using both the standard single¬ 
source approach with a modest level of smearing and our 
variational approach. 

The magnetic form factor, shown in Fig. HKb), is the 
most striking example. Here we see a clear difference in 
the quality of the plateau. For the correlation matrix 
approach, single-state dominance follows the current at 
tc = 21 immediately, allowing fits as early as ts = 22. 

For the standard approach, the excited states act to 
suppress the value of Gm at earlier time slices forcing one 
to wait until at least ts = 24 before an adequate x^of is 
obtained. However, the central values show a systematic 
upward trend following this time slice and consequently 
there is no indication that a plateau has been obtained. 

Similar conclusions can be drawn for the quadrupole 
form factor in Fig.|3](c). Here we find that both the corre¬ 
lation matrix and the standard approach give consistent 
values immediately following the current, but diverge as 
we move out to later time-slices. Again there is a clear 
systematic drift in the results obtained using the stan¬ 
dard approach and it would be difficult to select a fit 
region where the form factor can be determined with con¬ 
fidence. Given that we seek a region where the extracted 
form factor is constant over successive time slices, we can 
clearly see the improvement offered by the correlation 
matrix approach. 

In the case of the charge form factor, Fig. [UJa), the 
two methods are in closer agreement and display a sim¬ 
ilar quality in plateau. In either case, some Euclidean 
time evolution is required before a plateau is observed, 
but we find that the correlation matrix approach gives 
a systematically lower value following the current and 
plateaus a couple of time-slices earlier than the standard 
approach. 

Though the examples presented are selected to high¬ 
light the improvement using the variational method, we 
see significant improvement across all masses and form 
factors considered. Through this method we consistently 
find that one is able to obtain single-state dominance ear¬ 
lier and often with signihcant improvement in the qual¬ 
ity of the plateau. Through the use of the variational 
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(a) Charge form factor, Gc, for m,r = 296 MeV. 



(b) Magnetic form factor, Gm, for = 156 MeV. 



(c) Quadrupole form factor, Gq, for m,r = 296 MeV. 


FIG. 4. A comparison of p-meson form factors as a function of 
Euclidean sink time for a single level of smearing and our vari¬ 
ational approach. The data sets obtained from ratios of three- 
and two-point functions are offset for clarity. The (blue) cir¬ 
cles denote the results from the variational approach while 
the (red) squares illustrate traditional results using the stan¬ 
dard single-source method with a moderate level of smearing 
(35 sweeps). The vertical dashed line indicates the position 
of the current insertion. The fitted value from the variational 
approach has been included (shaded band) to highlight where 
the single source approach is consistent with our improved 
method. 


approach we are able to: 

• Isolate an eigenstate at earlier Euclidean times, 

• Insert the conserved vector current at earlier Eu¬ 
clidean times, 

• Fit the correlation function at earlier Euclidean 
times, 

• Observe robust plateau behaviour, 

• Identify large Euclidean-time fit windows, 

• Determine the form factors with significantly re¬ 
duced systematic errors, and 

• Determine the form factors with signihcantly re¬ 
duced statistical uncertainties due to the admission 
of analysis at earlier Euclidean times. 

C. Form Factors and Static Quantities 

In extracting form factors across a wide range of quark 
masses, it is important to note that for each mass there 
will be a slight change in the value of due to the 
variation in the mass of the hadron. To ensure that 
the comparison between quark masses and eigenstates is 
meaningful, we make use of the monopole ansatz to shift 
the extracted values for our form factors to a common 

. Fig. [5] demonstrates this shift for the pion. For the 
pion system, we shift values to a common value of 

= 0.1 GeV^, while for the p-meson system we select 

= 0.16 GeV^. These values are selected to minimise 
the shift for the extractions at the lightest quark mass. 
We use different values to ensure that we minimise the 
shift for each system. As the pion is signihcantly lighter 
than the p meson, a smaller value for arises naturally. 

Before we examine our form factor results, we note 
that all quantities presented are the quark sector contri¬ 
butions for unit-charge quarks. Here we choose to label 
these as the quark sector contributions to the positive- 
charge eigenstate of the corresponding iso-triplet. That 
is, the quark contribution is labelled as the u-quark sec¬ 
tor while the anti-quark contribution is labelled as the 
d-quark sector. As we are working with exact isospin 
symmetry, these quark sector contributions are equiva¬ 
lent and so we choose to present the u-quark sector only. 
As we are working in full QCD, the hadron form factor 
will have contributions from the sea quarks. However 
these contributions have been neglected in our calcula¬ 
tion. Thus to make meaningful comparison with exper¬ 
imental data, one should consider the iso-vector quanti¬ 
ties, 

(1f) = ^((1+)-(1-)), 

where (1±) labels the positive and negative members of 
the iso-triplet of the pion and p meson. However, un¬ 
der exact isospin symmetry one finds that this iso-vector 
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FIG. 5. An example of the shift applied to the pion form 
factor to ensure that all eigenstates across all quark masses 
are at a common . For the pion and its excitation, this 
shift is to = 0.10 GeV^, while for the rho meson and its 
excitation we shift to = 0.16 GeV^. 


FIG. 7. The unit-charge quark-sector contributions to the 
charge form factor, Gc, for the rho meson (blue circles) and 
its hrst excitation (red squares) at the common value = 
0.16 GeV^ 
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FIG. 6. The unit-charge quark-sector contributions to the 
charge form factor, Gc, for the pion (blue circles) and its 
first excitation (red squares) at the common value = 
0.10 GeV^. The dashed line represents the physical point. 


quantity defined above is exactly that given by the con¬ 
nected quark sector contribution. 


1. Charge Form Factors 


TABLE IV. The quark sector contributions to the charge form 
factor, Gc, for the ground state and excited state pi meson, 
at the common — 0.10 GeV^. Results are for unit charge 
quarks. 


ml (GeV^) 

'Utj- 

'Utt* 

0.3876(11) 

0.898(1) 

0.868(11) 

0.2647(9) 

0.884(2) 

0.828(16) 

0.1509(7) 

0.880(3) 

0.731(25) 

0.0811(6) 

0.846(5) 

- 

0.0260(10) 

0.870(16) 

- 


noted that this is consistent with quark model expec¬ 
tations where a hyperfine interaction of the form 
is repulsive when spins are aligned and attractive when 
spins are anti-aligned. In Fig. 10 we also include the ex¬ 
perimental value for the pion radius Q , which compares 
well with our determination. 

For the heaviest quark masses where well determined 
values are available for the excited states, it appears that 
a similar trend holds between our p* and tt* mesons. The 
relatively large radii for the excited states at light quark 
masses are interesting. Results are tabulated in Table 

M 


In Figs. [Hand[7] we display the charge form factor Gc 
for the TT and tt* mesons and p and p* mesons respectively. 
In both channels we observe a decrease in the charge form 
factor for the excitation which translates to the excited 
states having a larger spatial extent. Values are reported 
in Tables |TV] and |3 

To give us insight into the relative size of these states, 
we consider the charge radii, shown in Fig. [5] As was 
found in Q , the ground state vector meson is consistently 
larger than the corresponding pseudoscalar meson. It was 


2. Magnetic Form Factors 

The magnetic form factors for our p and p* mesons 
are illustrated in Fig. |H For both the ground and ex¬ 
cited state we observe very little variation in the value 
as rriT^ varies. As summarised in Table I VIII we clearly 
observe a significantly smaller value of Gm for the exci¬ 
tation. Though we would expect a decrease consistent 
with the decrease in the charge form factor for this state, 
the degree of suppression suggests that the magnetic mo- 
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TABLE V. The quark sector contribution to the charge form 
factor, Gc, of the ground state and excited state rho meson, 
at the common = 0.16 GeV^. Results are for unit charge 
quarks. 


ml (GeV^) 

Up 

Up- 

0.3876(11) 

0.801(3) 

0.765(12) 

0.2647(9) 

0.783(3) 

0.727(16) 

0.1509(7) 

0.772(3) 

0.749(14) 

0.0811(6) 

0.733(10) 

0.639(34) 

0.0260(10) 

0.685(22) 

0.640(58) 
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FIG. 8. Mean-square charge radii, (r^), for the positive-charge 
states of the tt (blue circles), tt* (green triangles), p (red 
squares) and p* (purple diamonds) mesons. The grey dashed 
line represents the physical point, while the light blue data 
point is the experimental value for the pion from PDG [j. 
The TT* and p* values have been offset for clarity. 


ment for this state is smaller than the ground state. 

In Table lVnil we list magnetic moments for these states 
where we invoke common scaling between the charge and 
magnetic form factors. We find a smaller moment for the 
p*. Our observation is consistent with results for meson 
magnetic moments using the relativistic Hamiltonian 
approach. In particular, their value for ~ 0.7 com¬ 
pares well with our result of 0.74(9) for the second lightest 
mass. The quark-mass flow of these magnetic moments 
is illustrated in Fig. [TOl 

In Fig. [TT]we show the g-factor for the p meson, pro¬ 
vided by the magnetic moment of the p meson in natural 
magnetons. Constituent quark model expectations sug¬ 
gest for a pure s-wave state, that Pp ~ 2. Our result of 
Pp = 2.21(8), taken from our lightest mass point, is larger 
than this expectation and suggests a non-trivial value for 
the quadrupole moment of the p associated with U-wave 
mixing. We observe a mild downwards trend of the p- 
factor with increasing quark mass suggesting that our 
results are compatible with the quark-model expectation 
of = 2 in the large quark mass limit. 

Our results agree in value and behaviour with the pre¬ 


TABLE VI. Mean-square charge radii ((r^)) for the positive- 
charge states of the tt and p mesons and their first excitations 
in units of fm^. 


ml (GeV^) 

TT 

P 

TT* 

P* 

0.3876(11) 

0.267(4) 

0.363(6) 

0.35(3) 

0.45(3) 

0.2647(9) 

0.306(7) 

0.405(8) 

0.48(5) 

0.55(4) 

0.1509(7) 

0.320(10) 

0.430(8) 

0.86(11) 

0.49(4) 

0.0811(6) 

0.425(16) 

0.531(27) 

- 

0.83(12) 

0.0260(10) 

0.349(51) 

0.670(68) 

- 

0.82(21) 


2.5 

2.0 


0.5 



FIG. 9. The unit-charge quark-sector contributions to the 
magnetic form factor, Gm, of the p"*" and p*^ mesons, at the 
common — 0.16 GeV^. 


vious quenched determination [5| an d dynamical study 
using background field methods |43| . In Fig. [TT]we in¬ 
clude the experimental determination of [l|. Within the 
literature, the majority of model calculations [dsl - l^ give 
a value of pp between 2.0 and 2.4 similar to our determi¬ 
nation of 2.21(8). 


3. Quadrupole Form Factors 

The quadrupole form factor for the p meson is shown 
in Fig. [T21 For the excited p* meson, the signal was too 
poor to extract a result. As was found in the quenched 
study of Ref. @ and the preliminary study of Ref. Q , the 
quadrupole form factor is negative in value. The value 
of Gq varies mildly in the heavy quark regime, however 
we observe a significant increase in the magnitude as we 
move to the lightest mass. Numerical results are provided 
in Table HXl 

In Fig. [13] we illustrate the quark mass dependence of 
the quadrupole moment. We also include the quadrupole 
moment extracted from the quenched data in Ref. @ . 
At heavier masses we find consistent values and find a 
slight increase in the magnitude in the direction of de¬ 
creasing quark mass. However, we see a rapid increase in 
the quadrupole moment magnitude at our lightest mass. 
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FIG. 10. The quark-mass dependence of the p"*" and p*+ mag¬ 
netic moments in units of the nuclear magneton. 



FIG. 11. The p-factor of the p meson is provided by the 
magnetic moment of the p meson in natural magnetons. The 
dashed line highlights the physical point. The li^t blue data 
point is the experimental determination of Ref. [l|. 


The value nearly doubles in comparison with the next 
lightest mass. This indicates the importance of light- 
quark dynamics to the underlying structure of the rho 
meson giving rise to significant contributions from the 
pion cloud. The dramatic variation observed warrants 
further investigation into the chiral dynamics of this 
quantity, especially into the role that finite volume ef¬ 
fects may have on the determination. 

In Ref. [1^ , through considerations of the most general 
free Lagrangian for a charged spin-1 system with mini¬ 
mal electromagnetic coupling, it was shown that there 
exists an explicit degree of freedom in the Lagrangian 
which can be parametrised by the p-factor. Conse¬ 
quently one finds that at tree-level the quadrupole mo¬ 
ment Qp = {1 — gpYjml. With our value of Pp = 2.21(8), 
the tree-level value for Qp ~ —1.21(8)e/mp. In Table DCl 
we report the quadrupole moment. To make contact with 
Ref. |53l ]. we report results in terms of the natural units 
of e/m^ where rup is the mass of the p meson observed 


TABLE VII. Unit-charge quark-sector contributions to the 
magnetic form factor of p and p* mesons, at the common 
= 0.16 GeV^. The natural magneton has been converted 
to constant units of nuclear magnetons pjv. 

ml (GeV^) 

Up (pjv) 

Up* (pn) 

0.3876(11) 

1.713(12) 

1.02(6) 

0.2647(9) 

1.779(21) 

1.41(12) 

0.1509(7) 

1.864(17) 

1.09(9) 

0.0811(6) 

1.978(56) 

1.27(14) 

0.0260(10) 

1.791(37) 

- 

TABLE VIII. Magnetic 

moments of the positive-charge states 

of the p and p* mesons 

in units of nuclear magnetons pv. 

ml (GeV^) 

Pp (pn) 

Pp* (pn) 

0.3876(11) 

2.138(15) 

1.33(8) 

0.2647(9) 

2.272(26) 

1.94(17) 

0.1509(7) 

2.414(23) 

1.45(12) 

0.0811(6) 

2.698(74) 

2.00(26) 

0.0260(10) 

2.613(97) 

- 



FIG. 12. The unit-charge quark-sector contribution to the 
quadrupole form factor, Gq, of the p meson at the common 
= 0.16 GeV^ 


on the lattice at each quark mass. At the lightest quark 
mass considered we find Qp = —0.733(99)®/mp indicating 
important contributions beyond tree level, driven by the 
fundamental strong interactions of QCD. 

From the perspective of a non-relativistic quark model, 
the quadrupole moment arises from an admixture of S 
and D-wave components in the wave function. Thus our 
non-zero quadrupole moment, even at heavy masses, in¬ 
dicates an important D-wave component to the p meson. 

In a next-generation calculation where excited-state 
signals are sufficiently accurate, the quadrupole moment 
can be used to determine the dominant contributions to 
the p(1450) and p(1700) wave functions. Thus, it of¬ 
fers the ideal tag to track these eigenstates with varying 
quark mass. In this way, one can determine if there is a 
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TABLE IX. Unit-charge quark-sector contributions to the 
quadrupole form factor of the p meson, at the common 
= 0.16 GeV^ 


ml (GeV^: 

) 

Up (fm^) 

0.3876(11) 


-0.0117(6) 

0.2647(9) 


-0.0134(9) 

0.1509(7) 


-0.0125(13) 

0.0811(6) 


-0.0161(22) 

0.0260(10) 


-0.0310(45) 

TABLE X. 

The quadrupole moment of the p"*" i 

meson in nat- 

ural units of elmi, where rrip is the mass of the p meson 

observed on the lattice at each quark mass, and 
of fm^ 

in hxed units 

ml (GeV^; 

) gr* 

Qp (fm^) 

0.3876(11) 

-0.362(20) 

-0.0146(8) 

0.2647(9) 

-0.368(25) 

-0.0171(12) 

0.1509(7) 

-0.313(32) 

-0.0162(17) 

0.0811(6) 

-0.392(53) 

-0.0219(30) 

0.0260(10) 

-0.733(99) 

-0.0452(61) 


reordering of these states in the light quark regime. 


VI. SUMMARY AND CONCLUSIONS 

We have established a general framework for the use 
of the variational approach in the evaluation of hadronic 
form factors. This approach can be used to system¬ 
atically eliminate excited state contributions to ground 
state matrix elements. It also allows one to evaluate the 
same quantities for excited states with no additional ef¬ 
fort. 

As we found in Ref. [l^, the use of optimised inter¬ 
polators obtained from the variational approach results 
in rapid isolation of the eigenstate, enabling earlier in¬ 
sertion of the probing current. Optimised interpolators 
at the sink result in rapid onset of robust plateau be¬ 
haviour enabling early and large Euclidean-time fit win¬ 
dows. Together these features act to reduce systematic 
errors through the suppression of excited state contami¬ 
nations and reduce statistical uncertainties through the 
ability to insert the current and establish fit windows ear¬ 
lier in Euclidean time. This approach, coupled with the 
large lattice volume and light quark masses, has resulted 
in an accurate determination of the pi and rho electro¬ 
magnetic form factors at the low . 

Our light quark-mass determination of the rho-meson 
g-factor, pp = 2.21(8), compares well with the experimen¬ 
tal result of [I| , but with significantly smaller uncertainty. 
This value is consistent with earlier lattice and model 
evaluations, which collectively prefer a g-factor slightly 
larger than the simple quark model estimate of 2. 



FIG. 13. The quark-mass dependence of the p-meson 
quadrupole moment. The blue circles are from the current 
analysis, while the red squares are the quenched results from 
Ref. @|. We see very similar behaviour between the results 
at the heavier masses. However, at our lightest mass we ob¬ 
serve a significant increase in the magnitude of the quadrupole 
moment indicating a signihcant role of the pion cloud in the 
underlying structure of the rho meson. 

As was found in the quenched calculation of Ref. @ , we 
obtain a negative value for the quadrupole form factor. 
The onset of significant chiral nonanalytic behaviour in 
the light quark-mass regime is also observed. 

Einally we have for the hrst time measured the elec¬ 
tromagnetic form factors for a light meson excitation. 
We Hnd that the charge form factors for these states are 
smaller than their ground state counterparts, consistent 
with expectations that these states should be larger in 
size. Eor the excited rho meson, we observed a signifi¬ 
cantly smaller value for the magnetic form factor and a 
smaller magnetic moment for the p*. Our observation 
of p.p*Jup = 0.74(9) supports the model prediction of 
Ref. gl]. 

Future work will investigate the hnite-volume correc¬ 
tions to these matrix element calculations drawing on the 
isolation of spin polarizations in the lattice simulations 
and effective field theory in relating the finite-volume ma¬ 
trix elements to those realised in Nature. 

For the generalised eigenvalue problem with explicit 
spin degrees of freedom on a finite-volume lattice, 
one can begin by quantifying the volume-induced non¬ 
degeneracies in the spectrum. For boosted eigenstates 
with spin, there is interplay between the spin polariza¬ 
tion direction and the momentum direction. For a fixed 
momentum direction, the reduced rotational invariance 
induced by the finite volume of the lattice gives rise to dif¬ 
ferent multi-particle dressings for different polarizations. 
This gives rise to a subtle variation in the energy of the 
eigenstates measured using different polarizations. To 
correctly account for this, a formalism to take the infi¬ 
nite volume limit of each spin-momentum combination 
considered is to be developed prior to combining the var¬ 
ious components to extract the form factors. 
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Two-point function analyses are in the exploratory 
phase [5^ and the effect is small. Thus, an exploration 
of the subtle nature of finite-volume corrections await 
a next-generation simulation. There the role of discon¬ 
nected quark loop contributions to the form factors can 
be ascertained. 


ACKNOWLEDGEMENTS 

We thank PACS-CS Collaboration for making their 
2 -1-1 flavor configurations available and acknowledge the 
important ongoing support of the ILDG. This research 
was undertaken with the assistance of resources at the 
NCI National Facility in Canberra, Australia, and the 
iVEC facilities at Murdoch University (iVEC@Murdoch) 
and the University of Western Australia (iVEC@UWA). 
These resources were provided through the National 
Computational Merit Allocation Scheme, supported by 
the Australian Government. We also acknowledge eRe- 
search SA for their support of our supercomputers. This 
research is supported by the Australian Research Coun¬ 
cil. 


[1] D. G. Gudino and G. T. Sanchez, (2013), 

arXiv:1305.6345 [hep-ph] 

[2] G. Lopez Gastro and A. Mar¬ 
iano, Phys.Lett. B517, 339 (2001) 

arXiv:nucl-th/0006031 [nucl-th] 

[3] M. Kotulla, J. Ahrens, J. Annand, R. Beck, 

G. Caselotti, et al., Phys.Rev.Lett. 89, 272001 (2002) 
arXiv:nucl-ex/0210040 [nucl-ex] 

[4] J. Beringer et al. (Particle Data Group), 

Phys.Rev. D86, 010001 (2012) 

[5] J. Hedditch, W. Kamleh, B. Lasscock, D. Leinwe- 
ber, A. Williams, et al, Phys.Rev. D75, 094504 (2007) 
arXiv:hep-lat/0703014 [HEP-LAT] 

[6] M. Gurtler et al. (QCDSF Collaboration), PoS LAT- 
TICE2008, 051 (2008). 

[7] C. Alexandrou, P. de Forcrand, and 

A. Tsapalis, Phys.Rev. D66, 094503 (2002) 
arXiv:hep-lat/0206026 [hep-lat] 

[8] S. Capitani, M. Della Morte, G. von Hippel, B. Jager, 

A. Juttner, et al., Phys.Rev. D86, 074502 (2012) 
arXiv: 1205.0180 [hep-lat] 

[9] S. Capitani, M. Della Morte, G. von Hippel, B. Jager, 

B. Knippschild, et al., PoS LATTICE2012, 177 (2012), 
arXiv:1211.1282 [hep-lat] 

[10] J. Bulava, M. Donnellan, and R. Sommer, 
JHEP 1201, 140 (2012), arXiv:1108.3774 [hep-lat] 

[11] T. Maurer, T. Burch, L. Y. Glozman, C. Lang, 
D. Mohler, et al, (2012), arXiv: 1202.2834 ]hep-lat] 

[12] B. J. Owen, J. Dragos, W. Kamleh, D. B. Leinwe- 
ber, M. S. Mahbub, et al., Phys.Lett. B723, 217 (2013) 
arXiv:1212.4668 [hep-lat] 

[13] B. J. Menadue, W. Kamleh, D. B. Leinweber, M. S. Mah¬ 
bub, and B. J. Owen, PoS LATTICE2012, 178 (2012). 

[14] B. J. Menadue, W. Kamleh, D. B. Leinweber, M. S. Mah¬ 
bub, and B. J. Owen, (2013), arXiv:1311.5026 [hep-lat] 

[15] B. Owen, W. Kamleh, D. B. Leinweber, S. Mahbub, and 
B. Menadue, PoS LATTICE2012, 173 (2012). 

[16] B. J. Owen, W. Kamleh, D. B. Leinweber, M. S. Mahbub, 
and B. J. Menadue, PoS LATTICE2013, 277 (2013), 
arXiv:1312.0291 [hep-lat] 

[17] S. Dinter, C. Alexandrou, M. Constantinou, V. Drach, 
K. Jansen, et al., Phys.Lett. B704, 89 (2011) 
arXiv: 1108.1076 [hep-lat] 


[18] H.-W. Lin, PoS LATTICE2012, 013 (2012), 

arXiv: 1212.6849 [hep-lat] 

[19] C. Alexandrou, PoS LATTICE2010, 001 (2010), 
arXiv:1011.3660 [hep-lat] 

[20] C. Michael, Nucl.Phys. B259, 58 (1985) 

[21] M. Luscher and U. Wolff, Nucl.Phys. B339, 222 (1990) 

[22] We note that an analysis accounting for finite-volume sys¬ 
tematic errors will need to first construct the spin eigen¬ 
states due to finite-volume induced non-degeneracies in 
the spectrum [s^ . 

[23] S. Boinepalli, D. Leinweber, A. Williams, J. Zan- 
otti, and J. Zhang, Phys.Rev. D74, 093005 (2006) 
arXiv:hep-lat/0604022 [hep-lat] 

[24] S. Boinepalli, D. Leinweber, P. Moran, A. Williams, 
J. Zanotti, et al, Phys.Rev. D80, 054505 (2009) 
arXiv:0902.4046 [hep-lat] 

[25] T. Draper, R. Woloshyn, W. Wilcox, and K.-F. Liu, 
Nucl.Phys.Proc.Suppl. 9, 175 (1989) 

[26] K. J. Kim and Y.-S. Tsai, Phys.Rev. D7, 3710 (1973) 

[27] R. Arnold, C. E. Carlson, and F. Gross, 

Phys.Rev. C23, 363 (1981) 

[28] S. J. Brodsky and J. R. Hiller, 
Phys.Rev. D46, 2141 (1992) 

[29] D. B. Leinweber, R. Woloshyn, and T. Draper, 

Phys.Rev. D43, 1659 (1991) 

[30] D. B. Leinweber, T. Draper, and 

R. Woloshyn, Phys.Rev. D46, 3067 (1992) 

arXiv:hep-lat/9208025 ]hep-lat] 

[31] S. Aoki et al., Phys.Rev. D79, 034503 (2009) 
arXiv:0807.1661 [hep-lat] 

[32] M. G. Beckett et al, 

Comput.Phys.Commun. 182, 1208 (2011), 
arXiv:0910.1692 [hep-lat] 

[33] G. Martinelli, C. T. Sachrajda, and A. Vladikas, 
Nucl.Phys. B358, 212 (1991) 

[34] S. Gusken, Nucl. Phys. Proc. Suppl. 17, 361 (1990) 

[35] M. S. Mahbub, W. Kamleh, D. B. Leinweber, 

P. J. Moran, and A. G. Williams (CSSM Lat¬ 

tice collaboration), Phys.Lett. B707, 389 (2012) 
arXiv:1011.5724 [hep-lat] 

[36] D. S. Roberts, W. Kamleh, and D. B. 

Leinweber, Phys.Rev. D89, 074501 (2014) 

arXiv:1311.6626 [hep-lat] 




16 


[37] B. J. Menadue, W. Kamleh, D. B. Leinweber, 

and M. S. Mahbub, Phys.Rev.Lett. 108, 112001 (2012) 
arXiv: 1109.6716 [hep-lat] 

[38] M. S. Mahbub, W. Kamleh, D. B. Leinweber, 

P. J. Moran, and A. G. Williams (CSSM Lat¬ 
tice Collaboration), Phys.Rev. D87, 011501 (2013) 

arXiv: 1209.0240 [hep-lat] 

[39] M. S. Mahbub, W. Kamleh, D. B. 

Leinweber, P. J. Moran, and A. G. 

Williams, Phys.Rev. D87, 094506 (2013) 

arXiv:1302.2987 [hep-lat] 

[40] J. J. Dudek, R. G. Edwards, M. J. 

Peardon, D. G. Richards, and C. E. 

Thomas, Phys.Rev. D82, 034508 (2010) 

arXiv:1004.4930 [hep-ph] 

[41] C. Gattringer, L. Y. Glozman, C. Lang, D. Mohler, 

and S. Prelovsek, Phys.Rev. D78, 034501 (2008) 

arXiv:0802.2020 [hep-lat] 

[42] G. P. Engel, C. Lang, M. Limmer, D. Mohler, 

and A. Schafer (BGR [Bern-Graz-Regensburg] 
Collaboration), Phys.Rev. D82, 034505 (2010) 

arXiv:1005.1748 [hep-lat] 

[43] A. Badalian and Y. Simonov, 

Phys.Rev. D87, 074012 (2013), 

arXiv:1211.4349 [hep-ph] 

[44] F. X. Lee, S. Moerschbacher, and 

W. Wilcox, Phys.Rev. D78, 094502 (2008) 

arXiv:0807.4150 [hep-lat] 


[45] J. de Melo and T. Fred- 

erico, Phys.Rev. C55, 2043 (1997) 

arXiv:nucl-th/9706032 ]nucl-th] 

[46] M. Bhagwat and P. Maris, 

Phys.Rev. C77, 025203 (2008), 
arXiv:nucl-th/0612069 [nucl-th] 

[47] T. Aliev, A. Ozpineci, and M. Savci, 

Phys.Lett. B678, 470 (2009), arXiv:0902.4627 [hep-ph] 

[48] H. Roberts, A. Bashir, L. Gutierrez- 

Guerrero, C. Roberts, and D. Wil¬ 
son, Phys.Rev. C83, 065206 (2011) 

arXiv: 1102.4376 [nucl-th] 

[49] M. Pitschmann, C.-Y. Seng, M. J. Ramsey- 

Musolf, C. D. Roberts, S. M. Schmidt, 
et al., Phys.Rev. C87, 015205 (2013) 

arXiv: 1209.4352 [nucl-th] 

[50] D. Djukanovic, E. Epelbaum, J. Gegelia, and 
U.-G. Meissner, Phys.Lett. B730, 115 (2014) 
arXiv: 1309.3991 [hep-ph] 

[51] C. S. Mello, A. Nunes da Silva, J. de Melo, and T. Fred- 
erico, Nucl.Phys.Proc.Suppl. 251-252, 62 (2014) 

[52] We note that an error was made in Ref. in converting 
the lattice form factors from natural units to fixed units. 

[53] T. D. Lee, Phys. Rev. 140, B967 (1965) 

[54] W. Kamleh, D. B. Leinweber, B. J. Owen, J. Wu, and 
R. D. Young, In preparation. 


